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Abstract 

The totally asymmetric simple exclusion process (TASEP) on the 
one-dimensional lattice with the Bernoulli p measure as initial con- 
ditions, < p < 1, is stationary in space and time. Let Nt(j) be 
the number of particles which have crossed the bond from j to j + 1 
during the time span [0, t}. For j = (1 — 2p)t + 2w(p(l — /o)) 1//3 i 2//3 we 
prove that the fluctuations of Nt(j) for large t are of order t 1 / 3 and we 
determine the limiting distribution function F w (s), which is a gener- 
alization of the GUE Tracy- Widom distribution. The family F w (s) of 
distribution functions have been obtained before by Baik and Rains in 
the context of the PNG model with boundary sources, which requires 
the asymptotics of a Riemann-Hilbert problem. In our work we arrive 
at F w (s) through the asymptotics of a Fredholm determinant. F w (s) 
is simply related to the scaling function for the space-time covariance 
of the stationary TASEP, equivalently to the asymptotic transition 
probability of a single second class particle. 



1 Scaling limit and main result 

The totally asymmetric simple exclusion process (TASEP) is, arguably, the 
simplest non-reversible interacting stochastic particle system. The occupa- 
tion variables of the TASEP are denoted by r/j, j e Z, rjj = means site j 
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is empty and rjj = 1 means site j is occupied. Since we plan to study the 
stationary space-time covariance (= two-point function), the particles move 
on the entire one- dimensional lattice Z. The stochastic updating rule is ex- 
tremely simple. Particles jump to the right and are allowed to do so only if 
their right neighboring site is empty. Jumps are independent of each other 
and are performed after an exponential waiting time with mean 1, which 
starts from the time instant when the right neighbor site is empty. 

More precisely, we denote by rj a particle configuration, rj G Q = {0, 1} Z . 
Let /: Q — > R be a function depending only on a finite number of r^'s. Then 
the backward generator of the TASEP is given by 

ljxv) = 5>a - - my (1.1) 

Here rf ,: > +1 denotes the configuration 77 with the occupations at sites j and 
j + 1 interchanged. The semigroup e Lt is well-defined as acting on bounded 
and continuous functions on Q. e Lt is the transition probability of the 

TASEP jug. 

Let \ip be the Bernoulli measure with density p, < p < 1, i.e., under p p 
the rj/s are independent and p p {r\j = 1) = P- From p.ljl it is easy to check 
that 

Pp(Lf) = (1.2) 

for all local functions /, which means that the Bernoulli measures are station- 
ary measures for the TASEP. In fact, these are the only translation invariant 
stationary measures [Hj. In the sequel we fix p, excluding the degenerate 
cases p = 0, p = 1, and start the TASEP with p p . The corresponding space- 
time stationary process is denoted by i]j(t), t G R, j G Z. Pta denotes the 
probability measure on paths t 1— >■ i](t) and Eta its expectation. The depen- 
dence on p is always understood implicitly. Note that the average current 
for the stationary TASEP is j(p) = p(l — p). 

As for any other stationary stochastic field theory the most basic quantity 
is the two-point function, which for the TASEP is defined through 

E T AM*)r/o(0))-p 2 = S(j,t). (1.3) 

For fixed t, S(j,t) decays exponentially in j. One has the sum rules 

s(j, *) = £ {E TA (vj(t)m(o)) - p 2 ) = p(l -p) = x(p), (1.4) 

-^^(j,t)=/(p)t = (l-2p)t. (1.5) 
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X 1 S(j,t) can be viewed as the probability for a second class particle to be 
at site j at time t given it was at j — initially, see e.g. [TH]. Thus 

S(j,t)>0, (1.6) 

which would not hold on general grounds. 
The next finer information is the variance 

^tf = X ~ l Y,3 2 S{3,t)-{{l-2p)tf. (1.7) 

Naively, one might expect that a(t) = y/i, arguing that the second class 
particle moves random walk like. As noticed in , in a purely perturbative 
argument, a(t) is likely to grow faster than \fi. The proper scaling form was 
firmly established in [3] with the result 

a(t) - a oX 1/3 t 2/3 (1.8) 

for large t. x 1 ^ 3 follows on dimensional grounds, while the prefactor ao has to 
be determined numerically. In fact, ao = 2.0209 . . . which is a consequence 
of the result reported here together with |21j . 

Forster, Nelson and Stephen [H] consider, as a particular case of the fluc- 
tuating Navier-Stokes equation, the stochastic Burgers equation 

5 B = s(-" , + "s" +{ ) (L9) 

with v > and £ space-time white noise, which is a sort of continuum 
stochastic partial differential equation version of the TASEP. They obtain 
the dynamical exponent z = 3/2 which corresponds to the 2/3 of (jl. 8)1 . Kar- 
dar, Parisi and Zhang [K] study surface growth which for a one- dimensional 
substrate reduces to ()1.9|) with u being the gradient of the height function. 
By more refined arguments they confirm z = 3/2 in one space dimension. 
Since then many approximate theories have appeared, see e.g. (THl HB| for 
a more complete discussion. The only one which survives the test is the 
mode-coupling theory, which is a nonlinear equation for S(j,t) [3]. A care- 
ful, rather recent, numerical study |3] of this equation yields surprisingly 
good agreement with the exact two-point function in the scaling limit [2*T] . 
The power law (jl.8|) strongly suggests the scaling form 

S(j,t) - \x(?X l/3 t 2/3 )- l 9 : c {{3 - (1 - 2p)t)(2 X 1 /3 t 2/3 r i) (L1Q) 

for large t and for j — (1 — 2p)t = 0(t 2 ^ 3 ) with the scaling function g'lJl 
independent of p. Our main result will be to prove a version of (jl.lOj) with 
a reasonably explicit expression for g" c . 



3 



The scaling function g sc appeared already in the context of the variance 
of the height differences in the polynuclear growth (PNG) model, where g sc 
is somewhat indirectly determined by a set of differential equations, which 
were discovered by Baik and Rains j2] , see Appendix \K\ These differential 
equations are solved numerically in [2*T] . where also a plot of g" c is displayed. 
Thus one important consequence of our main result is to establish that in the 
scaling limit the PNG model and the TASEP have the same scaling function 
for their covariance. Such a property is expected for a much larger class of 
one-dimensional driven lattice gases. For example, if instead of the TASEP 
we allow for partial asymmetry, to say a particle jumps with probability p to 
the right and 1 — p to the left, p ^ 1/2, then (jl.lU|) should still hold provided 
(1 — 2p)t is replaced by {2p — 1)(1 — 2p)t. The general formulation of the 
universality hypothesis for one-dimensional driven lattice gases is explained 
in [IB], see also [T9~] . Viewed in this context our main result asserts that the 
TASEP and the PNG model are in the same universality class. 

The issue of universality is certainly one strong motivation for our study. 
At first sight PNG and TASEP look very different, while when viewed prop- 
erly they are in fact not so far apart. The interpolating family of models is 
the TASEP with a discrete time updating rule. Its extreme limits are the 
PNG model on one side and the continuous time TASEP on the other side, 
sec [19 . Given the scaling limit for the PNG. it is not surprising to have the 
same result for the TASEP. However, it turns out that the method used in 
proving the analogue of (jl.lOj) for the PNG model does not generalize to the 
TASEP, which for us is an even more compelling reason to investigate the 
TASEP. At a certain stage the proof in [21] uses that the two-dimensional 
Poisson process is invariant under linear scale changes, which is a property 
special to the PNG model. For the TASEP we have to develop a novel 
method which will be rather different from [2] and uses non-intersecting 
line ensembles. In fact, our expression for g" c has an appearance quite unlike 
to the one discovered by Baik and Rains. It requires an argument that both 
expressions are in agreement, see Appendix 1X1 

To state our main result we have to first reformulate the TASEP as a 
growth process by introducing the height function h t (j) through 



t > 0, where N t counts the number of jumps from site to site 1 during the 
time-span [0,t]. Note that h t (j) — h (j) = 2N t (j), where N t (j) counts the 
number of particles which have crossed the bond from j to j + 1 during the 




for j > 1, 
for j = 0, 
for j < — 1, 



(1.11) 
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time span [0, t]. By stationarity one has 

Eta(Mj)) = 2p(l - p)t + (1 - 2p)j. (1.12) 

Since h t (j + 1) — h t (j) = —2r]j + i(t), the variance of the height must be 
simply related to S(j,t). 

Proposition 1.1. Let A be the discrete Laplacian, (Af)(j) = f(j + 1) + 
III 1) 2. /•(./). Then 

SS(j,t) = (AEta(M-) -E TA (h t (-))] 2 ))(j). (1.13) 

The proof can be found in Proposition 4.1 of [T§] . 
We introduce the family of distribution functions 

F w (s,t) = P TA ({(l-2 X )t + 2w(l-2p) X 1/3 t 2/3 -s X 2/3 t 1/3 

<h t ([(l-2p)t + 2w X 1/3 t^\)}), (1.14) 

where [^J denotes the integer part of x. Here the height is evaluated at 
(1— 2p)t, which is determined by the propagation of a tiny density fluctuation, 
plus a in comparison small off-set of order t 2//3 , while the distribution function 
is centered at ExA^tO')) w ^ n J = L(l — 2p)t + 2wx 1 ^ 3 t 2 ^ 3 \ and has an 
argument, — s, which lives on the scale x 2 ^ 3 t 1 / 3 . 

As to be shown, the distribution function F w (s,t) converges to a limit as 
t —>■ oo. The limit will be expressed in terms of a scaling function g and the 
GUE Tracy- Widom distribution function -Fgue(s). The latter can be written 
as a Fredholm determinant in L 2 (R), 

F G ue(s) = det(l - P K Ai>s P ) (1.15) 

with P the projector operator on [0, oo) and K Ais the integral operator with 
the Airy kernel shifted by s, i.e., 

K Ai , s (x,y) = / d\Ai(\ + x + s)Ai(\ + y + s). (1.16) 
Jr+ 

Define the functions 





fiR- 


dze wz K AitS {z,x)e ws 










</r,_ 


dze wz Ai(y + z + s) 








Ps(x,y) = 


(1- 


- P K Ai>s P Q )- l {x,y) 
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and the scaling function g by 



g(s,w) = e 3 



I„,3 



dxdye w{x+v) Ai(x + y + s) 



+ / dxdy^ WtS (x)p s (x,y)^ w>s (y) 



;i.i8) 



Our main theorem asserts the limit of the family of distribution functions 
F w (s,t). 

Theorem 1.2. Let -Fque an d g defined above. Then for fixed c\ < c<i one 
has 

lim / F w (s,t)ds = F GVE (c2 + w 2 )g(c 2 + w 2 ,w)-F Gl j E (c 1 +w 2 )g(c 1 + w 2 ,w) 
J Cl 

(1.19) 

pointwise. 

Corollary 1.3. The limiting height distribution function F w (s) is given by 

F w (s) = -^- s (F GlJE (s + w 2 )g(s + w 2 ,w)). (1.20) 

For the PNG model Baik and Rains obtain the limiting height distribution 
function denoted by H(s + w 2 ; w, —w) in Defintion 3 of [2 . It has the same 
structure as F w (s). Only the scaling function g is given as the solution of a 
set of differential equations, see Appendix 1X1 

The two-point function of the TASEP carries information on the variance 
of height differences, see (11.13)1 , while Theorem 11.21 provides the full family 
of distribution functions. In this sense ()1.19j) is a stronger result than (jl.K)j) . 
On the other hand, the limit (jl.lDj) for the distribution function asserts only 
the weak convergence of the corresponding probability measures, while from 
()1.13|) we infer that for the space-time covariance the convergence of second 
moments would be needed. If we assume a suitable tightness condition on 
F w (s, t), then 

hm J s 2 dF w (s,t) = J s 2 dF w (s) =g sc (w), (1.21) 
which together with Proposition 11.11 yields 

lim 2 X 1/3 t 2/3 ^(L(l " 2p)t + 2w X 1/3 t 2 / 3 \,t) = ~ X <&H (1.22) 
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when integrated against an arbitrary smooth function in w, in agreement 
with the claim (jl.lOj) . Tightness is also missing in the analysis of the PNG 
model. 

Over the recent years there has been a considerable interest in scaling 
limits for the TASEP. Slightly more general than here, one considers an initial 
measure which is Bernoulli p_ in the left half lattice Z_ and Bernoulli p + in 
the right half lattice 7L+. The initial step, p_ = 1 and p + = 0, is studied by 
Johansson [10J by mapping the TASEP to a last passage percolation problem. 
For general p + and p_ such a map is still possible and yields a last passage 
percolation problem with boundary conditions [15] . Through the Robinson- 
Schensted-Knuth (RSK) correspondence one then obtains a line ensemble 
with boundary sources. This line ensemble is determinantal, in fact only a 
rank one perturbation of the line ensemble with tie-down at both ends. We 
refer to [5J, where a similar construction is carried out for the line ensemble 
corresponding to the discrete time TASEP. There is also a link to the work 
by Baik, Ben Arous, and Peche pQ, who study rank r perturbations of the 
complex Gaussian sample covariance matrices. Viewed from this perspective 
the stationary TASEP is singular, which is partially overcome by the shift 
argument, see also [H]. But even then, in the resulting matrix elements there 
is still a delicate cancellation which tends to hide the asymptotics. The 
technique of line ensembles can be used also for the investigation of multi- 
point statistics [213 E], which however will not be needed in our context. 

In computer simulations mostly deterministic flat initial conditions are 
adopted, which translates to the initial particle configuration . . . 010101 . . . 
of the TASEP. As established by Sasamoto [22], see also [Zj, the single point 
statistics in the limit of large times is then given by the distribution of the 
largest eigenvalue of the GOE of random matrices and thus different from 
the distribution obtained in this contribution. For the PNG model the cor- 
responding result is proved prior by Baik and Rains [Sj, see also [Hj. 

Our paper is divided into two parts. The first part is a fixed t discussion 
of F w (s, t) with the goal to obtain a manageable expression. The second part 
is devoted to the asymptotic analysis. In the Appendices we establish that 
our expression for F w (s) agrees with the one of Baik and Rains, provide some 
background on the determinantal fields turning up, and explain how the line 
ensemble is related to the Laguerre kernel. 

Acknowledgements. 
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teraction and Growth in Complex Stochastic Systems". Michael Prahofer 
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2 Map to a directed polymer 

The statistics of the height function h t (j), restricted to the cone 
{jj^llil < h}, can be represented through a directed last passage percola- 
tion, see JH] • For this purpose, in the initial configuration, let ( + + 1 be the 
location of the first particle to the right of (and including) 1 and let — 
be the location of the first hole to the left of (and including) 0. Therefore 

C+ are independent and geometrically distributed, Q(C- — n) — (1 — p)p n , 
Q(C+ = n) = p(l — p) n , n — 0, 1, — In addition we define the family of in- 
dependent exponentially distributed random variables w(i,j), i,j > 0, such 
that w(i,j) has mean 1 for i,j > 1, w(i,0) has mean (1 — p)^ 1 for % > 1, 
w(0,j) has mean p -1 for j > 1, and w(0,0) = 0. The joint distribution of 
the random variables ( = (£+,£_) and {w(i, j),i, j > 0} is denoted by Q. 
These exponentially distributed random variables are linked to the TASEP 
in the following way: w(( + + £, 0) is the £-th waiting time of the first particle 
to the right of and w(0,(- + £) is the £-th waiting time of the first hole 
to the left of 0, £ = 1, 2, . . .. To describe the other w(i,j)'s we label in the 
initial configuration the particles from right to left such that the first particle 
to the right of has label 0. Then w(i,j), i,j > 1, is the j-th waiting time 
of particle i, where the first waiting time refers to the instant when the i-th 
particle is at lattice site — i + 1. 

For given ( let 



The w^(i, j)'s are used as local passage times in a directed last passage per- 
colation. Let us consider an up/right path u on N 2 with a finite number of 
steps. To it we assign the passage time 




forl<z<C+,j = 0, 

fori = 0,l< j<C-, 
otherwise. 



(2.1) 




(2.2) 
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Then the last passage time from point (0, 0) to point (m, n) is given by 

G(m,n) = max TiuS). (2.3) 

w:(0,0)-+(m,n) 

Here the maximum is over the up/right paths which start at (0,0) and end 
at (m, n). 

Proposition 2.1. [TH] With the above definitions 

Q({G(m, n) < t}) = Pta(W + n < h t (m - n)}). (2.4) 

G(m, n) can also be viewed as a growth process. We introduce the corre- 
sponding height function h(j,r), j € Z, r e N, through 



Mj,r) = 0, for |j| >r, (2.5) 

G((r - 1 +j)/2, (r - 1 - j)/2), if (-1)^ = -1, 
. G((r - 2 + j)/2, (r - 2 - j)/2), if (-1)^ = 1, 
for \j\ < t. 



h(j, r) 



By Proposition 12. 1[ 

U " (Pta({t -2 </*(;)}), if(-l) T+J = l. 1 ' 

/i(j, r) is not such a convenient quantity and we modify it by allowing an 
error of order 1. We display the (^-dependence of G(m,n) explicitly through 
G^(m, n). In particular, G°(m, n) is the random variable obtained by setting 

C+ = o = C- 

Proposition 2.2. Uniformly in the endpoint one has 

Q({|G c (m,n) - G°{m,n)\ > u}\() < C+e' u/{1 ' p) + (-e~ u/p . (2.7) 

Proof. We fix the endpoint (m, n). Let T^(u) be the passage time for w^(i,j) 
and let uj^^ be a maximizing path from (0,0) to (m,n). Then G^(m, 
T^Lx), G°(m,n) = T°« ax ). One has 



,n 



G°(m,n)-G<(m,n) = T°« ax ) - + T<« ax ) - T^Lx., 

< T°« ax )-T^« ax ) 

C+ C- 

< ^w(i,0) + ^^(O, j), (2.8) 
i=i j=i 
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where in the first inequality we used that a;^ ax is a maximizer of TV Similarly 

G°(m,n)-G<(m,n) = T°«J - T>LJ + T>LJ - T<(u4j 

> tVLx)-^L x )>o. (2.9) 

Combining Ij2"3|) yields 

C+ C- 

Q({|G c (m,n)-G°(m,n)| >n}|C) < <Q({$>(i, 0) + £ w(0, j) > u}\C 

i=i j=i 

< C+e~ u/(1 ~ p) + C-e Wp . (2.10) 

□ 

Definition 2.3. h°(j,r) is the height function as given in \2. where 
G(m, n) is replaced by G°(m, n) with the corresponding passage times Wo(i,j). 

It follows from the identity 1)2.6)1 and Proposition 12.21 that for e > 0, 
independent of t, one can choose d(e) such that 

(1 - e)Q({h°(j, r)<t- d(e)}) < F w (s, t) < Q({h°(j, t) <t+ d{e)}) + e, 

(2.11) 

uniformly in j, r. Therefore, setting 

*£M) = Q({h°([(l-2p)t + 2w X 1/3 t 2 / 3 \, (212) 

L(l - 2x)t + 2W(1 - 2^/3,2/3 - ^/3 t l/3j ) < t}) ? 

one has 

(l-e)F°(s,t-d(e)) < i^fot) < F°(s,t + d(e)) + e. (2.13) 

Since £ — > oo, we conclude that Theorem II .21 is implied by Theorem 12 . 41 st at ed 
below. 

Theorem 2.4. For fixed c\ < c<i the following limit holds 

lim f * F°( S ,£)ds = r F w ( S )ds. (2.14) 
The remainder of the paper deals with the proof of Theorem 12.41 
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3 The Laguerre line ensemble with boundary 
sources 

Let us consider the directed polymer in the general case of independent 
w(i,jYs with exponential distribution of mean 1/ay, i, j > 0. The directed 
polymer is determinantal (a notion which will be explained below) provided 
dij = o,i + bj > 0. For the case of Theorem 12 .41 one has to deal with Wo(i,j) 
and the obvious choice would be aj = \ — (1 — p)5i,o, bj = \ — p5j$. The 
corresponding directed polymer fails to be determinantal on two accounts: 

(i) ao + bo = whereas it should be striclty positive, 

(ii) formally w(0, 0) is uniformly distributed on R + , while in actual fact 
w (0,0) = 0. 

Our strategy is to first discuss the line ensemble for Oj = \ + (a — §)<^,o, 
&j = | + (b — !)<5j,o, a > 0, 6 > 0. This is the task of the current section. In 
the following section we will show that the case w(0, 0) = can be deduced 
from a shift argument. The resulting expressions will then be analytically 
continued to a — p — |, b — | — p. In this limit we recover w (i,j), which is 
required for Theorem 12.41 

Definition 3.1. Let w a ^(i, j) , i,j G N, be a family of independent exponen- 
tially distributed random variables such that 

H^a^jT 1 = 1 + (a - \)ko + ( b ~ |)^,o (3-1) 
with < a, b < \. 

With w a ^(i,j) as in Definition 13.11 let T(u) be as in ()2.2|) with w^(i,j) 
replaced by w a: b(i,j) and let 

G(m,n) = max TYw), (3.2) 

ui: (0,0)^ (m,n) 

compare with ()2.3|) . We define the height function h(j,r), j G Z, r G N, 
through ()2.5|) . It can also be generated by the following growth process, 

h(j,0) = 0, (3.3) 

{max{/i(j - 1, r), /i(j + 1, r)} 
+ u>«, 6 ((t + j)/2, (r - j)/2), if (-l) J+r = 1, 
h(j,r), if (-1)^ = -1, 

for \j\ < t + 1, 
/i(j,r + l) = 0, for |j|>r + l. 
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Figure 1: The growth dynamics associated with the TASEP directed last 
passage percolation. 



The dynamics is best visualized by extending h(j, r) to a function over 
R, through h(x,r) = h(j,r) for j — | < x < j + |, see Figured] Then 
alternately there is a stochastic and deterministic up-date. In the stochastic 
up-date mass is added to the current height h(x, r) according to w a ^(i,j), see 
(jH.Hj) . In the deterministic up-date down-steps move one unit to the right and 
up-steps one unit to the left. Thereby parts of the up-dated h may overlap. 
The maximum rule means that the excess mass in the overlap is annihilated. 

Underlying the growth process one may construct the corresponding 
Robinson-Schensted-Knuth (RSK) dynamics [H], which in our case simply 
means that the overlap annihilated in line I is copied to the lower lying line 
£ — 1. In formulas we set 

ht{j,0) = 0, (3.4) 

{ht-i(j,r) - hi(j,r) 
+ mm{h e (j - l,r),h e (j + l,r)}, if (-1)^' = 1, 
ht-i(j,r), if (-1)^ = -1, 

with the line label t = 0, — 1, . . .. 

The purpose of the RSK construction consists in having, for fixed r, a 
manageable statistics of the collection of points {he(j,r), £ G Z_, \j\ < t, 
ht{j,r) > 0}. To describe their statistics directly without recourse to the 
stochastic dynamics we first have to define admissible point configurations. 
Let {xj,j = —n, . . . , 0} be points on [0, oo) ordered as < x_ n < . . . < xo- 
We say that {xj,j = —n, . . . , 0} -< {x'j,j = —n, . . . , 0} if Xq < x' Q , Xj < x'j < 
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Figure 2: A non- intersecting line ensemble at r = 5. 



Xj+i for j = —n, . . . , —1. Admissible point configurations are then 

h e {±r,r) = 0, (3.5) 
{hiU,r),ieZ-} -< {h e (j + l,T),eeZ_}if |j|<rand(-l)^ = -l, 
{h e {j,r),£eZ_} y {h e (j + l,r),£ e Z_} if \j\< rand (-1)^ = 1. 

As with the growth dynamics, the order -< and y can be visualized by 
extending hi(j, r) to R, by setting hi(x, r) = ht(j, r) for j — | < x < j + |. 
Then f!3.5|) means that the lines he(x,r) do not intersect when considered as 
lines in the plane, see Figure 12 

To a given point configuration, alias line ensemble, one associates a 
weight. It is the product of the weights for each single jump. Let us use 
S as the generic symbol for a height difference. Then the up-step ho(—T,r) 
to ho(— t+1, r) has weight e -6 ' 5 ' and the down-step h (r — 1, r) to h (r, r) has 
weight e _a ' 5 '. All other jumps of the form hi(j, r) to h^(j + 1, r) have weight 
e ~l 5 l/2_ Note that the weights are assigned by reading the vector a from right 
to left and the vector b from left to right. The total weight is normalized to 
become a probability. This probability measure is called the Laguerre line 
ensemble with boundary values a, b. It agrees with the probability measure 
at growth time r obtained from the growth dynamics ()3.4j) together with the 
RSK construction. 

It is convenient to think of {ht(j, t),£ e ZL, j G Z} as a point process on 
Z x (0, oo), where j is referred to as time and hi as space. The corresponding 
random field is then 

0r(j,y) = ^5(^(j,r)-y), y>0. (3.6) 

According to our construction, at Zx {0}, i.e. at y — 0, there are an infinite 
number of points. However, the point measure refers only to points with a 
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strictly positive y coordinate. In fact, <f> T (j,y) is supported by 

L(2r-1)/4J 

E (2^-4j-l) (3.7) 

3=0 

points, ['J denoting the integer part. The point process <f) T (j,y) is determi- 
nantal, in the sense that is has determinantal moments. This means that 
there exists a kernel K a T ^(j,y;j',y') such that for a time-ordered sequence 
ji < • • • < im and arbitrary space-points y±, . . . ,y m > the m-th moment of 
(fi T is given by 

E (n^ fc '^)) = det ( K a,b(j*>yk'>jk'>yv))i<k,k><m- ( 3 - 8 ) 

fc=l 

For the two-point function of the TASEP we need the statistics of the 
random field <p T (j, y) only at fixed time j and in the remainder of this section 
we will provide an expression for K^(j,y;j,y'), y,y' > 0. There is no 
difficulty in principle to extend the construction and our results to unequal 
times. 

The distinction between odd and even j + r is slightly cumbersome and 
we restrict to odd r, even j by setting 

r = 2m+l, j = 2d. (3.9) 

In L 2 (R) we define P + as projection onto R + , P + + P_ = 1. We also 
introduce the operators T + , T_ with integral kernels 

T+(x,y) = e-(*-yy 2 e(x-y), 

T_(x,y) = e-^/ 2 Q(y-x), (3.10) 

where Q(x) — 1 for x > and Q(x) = for x < 0. In Fourier space T + is 
multiplication by (\ + ik)~ 1 and T_ by (| — ik)' 1 . Eigenfunctions of T + , T_ 
are the exponentials V'o(^) — e~ ax . For a < \ one has 

T + ^a(x) = -^—i> a {x) (3.11) 

2 - a 

and for a > —\ one has 

T.i) a (x) = — |— ^ a (x). (3.12) 
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For \d\ < m let, as an operator in L 2 (R), 

K m4 = LP_R (3.13) 

with 

Ij = r prn+drp—{m—d) jr, _ rp-(m+d)rpm-d ^ -j^-j 

It follows that, for a G (—1/2, 1/2), 

{R^- a ){x) = Z(a)^ a (x), (3.15) 
(L*^ a )(x) = Z(a)-^ a (x) 

where 

Z(a) = 2 + - . (3.16) 

(\ - a) m - d 

For later use we provide a representation of the kernel of K m ^. This kernel 
has a singular part, which is concentrated on the diagonal {x = y}. In the 
computations only the regular part will be used, hence only it is displayed. 
Since in Fourier space T + , resp. T_, is the operator of multiplication by 
(| + iky 1 , resp. (| — ik)^ 1 , with the change of variable \ — ik = z + p we 
obtain an integral expression for the regular part of L and i?. Let T p be a 
path around the pole p oriented anti-clockwise. Then the regular part of the 
kernels are 

^,y) = ^e {1/2 - p){x - y) irnAx-y), x-y>o, (3.17) 



where 

{z + p) m - d 



£\m+d ' 



(3.18) 



and similarly 



^^ie^'^/^fe-a;), x - y < 0, (3.19) 

where 

LAv -x) = l p ^ Z{y - X) { \/^1T - (3- 2 °) 

Let P u be the projection onto [u, oo). Then for any u > one has 

(P u K m>d P u ) (x, y) = G(x - u) [ dwL(x, w)R(w, y)Q{y - u). (3.21) 
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As explained in Appendix the regular part of K m d is a similarity 
transformed Laguerre kernel. Hence we refer to K m>d also as Laguerre kernel. 
It has the following properties, which are proved in Appendix El One could 
also arrive at an equivalent kernel by taking the exponential limit of the 
geometric case studied by Okounkov in [T7j . 

Proposition 3.2. Letu > 0. Then \\P u K m ^P u \\ < 1. In addition, for a > 0, 

Pu(l - K m>d )ip a G L 2 (R), P u (l - K m4 )*^ a e L 2 (R) (3.22) 
with a norm uniformly bounded in u. 

With these preparations we state the relation between the equal time 
kernel of (|3.8|) and the Laguerre kernel. 

Proposition 3.3. Let < a, b < |. Then for \d\ < m, x, y > 0, one has the 

identity 

K^ +1 \2d,x-2d,y) = K a m b d (x,y) (3.23) 

= K m4 (x,y) + -=—(1 - K m>d )iJ) b (x)(t - K mtd )*tfj a (y) 

with 

1 /l — 2a\ m (\ — 2b\ m (\ 9 \- d /l ,n\~ d 



4 Shift construction 

Let us consider the Laguerre line ensemble with boundary values a, b > 
and denote its weight by W a) &. Under W a ,b we want to study the weight of 
{ho(j, T ) < u } denoted by W a ,b({ho(j, r) < u}). More general events could 
be investigated, but there is no need in our context. We set w a ^(0,0) = v 
and recall that its weight is given by e~ v ( a+b \ v > 0. We display the explicit 
dependence of W a> b on v as W 0) &( , ,i>). 

From the construction of the Laguerre line ensemble one has, for v > 0, 
v + 5 > 0, the shift 

W a , 6 ({/ioO', r) < u}, v + 5) = e-^WatdhoU, r) + 5 < u}, v) (4.1) 

and differentiating 

^W 0j6 ({/io(i, r) < u}, v) = -(a + b)W a , b ({h (j, r) < u}, u) 

-H-Wa,b({h (j,T) <u},v). (4.2) 
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Since W a ,&(-) = / °° dvW ajb (-,v), by integrating in v, 

W Oj6 ({/i (i,r) <u},0) = ^W a , b ({h (j,r)<u}) 

+ (a + b)W a , b ({h (j,r) <«}). (4.3) 

Note that the left hand side is the weight for u> a ,&(0, 0) = 0. 

Let Z a)b (u) = W a ,b{{h {j,r) < oo},v) and Z af) = §™ dvZ a>b (v). Then, 
taking u — ► oo in (j4.3|) . 

^a,b(0) = (a + &)£«,&, (4.4) 

Z a , fe given in JS23). 

Let Pq' 6 be the probability for the Laguerre line ensemble in case 
w a ,&(0, 0) = and P Qj fe the one in case w a , b (0, 0) is exponentially distributed 
with mean (a + &) _1 , as in Definition IH.11 Then, by (|4.3|) and ()4.4j) . 

P a ' 6 ({/io(i,T)< W }) = ^(^^({^o(j,r)< M }) 

+(a + 6)P aj6 ({/ioO» <«})) (4-5) 

for u > 0. 

For determinantal processes probabilities as on the right hand side of 
()4.5|) are easily computed with the result 

Pa,b({h (j, r) < u}) = det(l - P u K a J d P u ), (4.6) 

where, as before, r = 2m+l, j = 2d, and P u projects onto the interval [u, oo). 
The determinant is regarded in L 2 (P) and the identity ()4.6|) makes sense only 
for u > 0. Thus we fix w > throughout. Since by (|3.23j) P u K°^ d P u is a rank 
one perturbation of P u K m ^P u and since 1 — P u K m ^P u is invertible, compare 
with Proposition I3.2| one arrives at 

det(l - P u K a m h d P u ) = det(l - P u K m4 P u )(a + b)G a > b (u) (4.7) 

with 

(a + b)G a > b (u) = l-^(^(t-K m4 )P u (l-P u K m4 P u y 1 P u (l-K m4 )^ b ) 

^a, b 

(4.8) 

with (•, •) denoting the inner product in L 2 (H). 
We also define 

F(u) = det(l - P u K m4 P u ) (4.9) 

and supply the m, d dependence of G a,b (u) and of .F(m) when needed. We 
summarize as 
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Proposition 4.1. Let w a ^{i, j) be as in Definition ^. 11 except for w aj b(0, 0) 
for which we set w a ^(0, 0) = 0. Let h be the corresponding top line as given 
in J23) ; (H3P- Then for u > 

P a '\{ho(j,r) < u}) = ^(F(u)G a >\u)) + F(u)(a + b)G a > b (u), (4.10) 

where F(u) is given in \4-9jj and G a,b (u) in fl^.&l ). 

5 Analytic continuation 

We have to extend the validity of (j4.10J) from < a, b < 1/2 to a + b = 0, 
which will be achieved by proving that both sides of ()4.1U|) are analytic. 

Proposition 5.1. The map (a,b) i— > Pq b {{ho{j, t) < m}) zs real analytic for 
a,b> -1/2. 

Proof. h (j,r) is measurable with respect to the cx-algebra generated by 
w(i,0), w(0,j), i,j = 1, . . . ,r. Let 

• • • ,Cr,^l, • ■ • ,^r) 

= PS* 6 ({/^o(j,r) < u}\w(i,0) = Zi,w(aJ) = Vj,i,3 = h ■■■,?) (5.1) 

as conditional probability. Clearly V u does not depend on a, b and < 14 < 1. 
Then 

P a ' b ({h (j,r)<u}) = [ f[(d^ + a)e^ 1+2 ^ 2 ) (5.2) 

T 

X II (<Ms + %~ (1+26)W2 )K(£l, • • • ,tr,Vl, ■ ■ -,Vr), 
k=l 

which by inspection is real analytic for a,b > —1/2. □ 

Proposition 5.2. Let u > and let G a ' b (u) be given by fl^.&| ). Then (a, b) \— > 
G a,b (u) extends to a real analytic function for a,b E (—1/2, 1/2). 

Proof. We repeat Eq. (l4~£j) . 

(a + b)Z a ^ b G a \u) = Z afi -{^ a ,(t-K m . d )P u (t-P u K m4 P u )- l P u {t-K m4 )^ h ) 

(5.3) 

with ip a (x) = e ax . 

First remark that (a + b)Z a ^ is real analytic for a, b G (—1/2, 1/2). On 
the other hand, in ()5.3j) the first (thus also the second) term diverges as 
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a + b — > 0. Thus we have to find another representation of G a,b such that 
both terms remain finite in the a + b —>■ limit. From each term we subtract 
the quantity (ip a , P u (t — K m4 )ip b ) and obtain 

r.h.S. Of Q = (Z a ,b- (lpa,PM) + Ma, PuK m , d ll> b ) (5.4) 

-(iP a , Q u (l - K m>d )P u (l - PuK^P^P^l - K m , d )^ b ) 

where Q u = 1 — P u . With this rearrangement one singles out the divergence 
which now are in Z ajb and (ijj a , P u ipb) only. We discuss the analytic continua- 
tion from a, b G (0, 1/2) to a, b G (—1/2, 1/2) for the three terms separately, 
where we use the properties (see proof of Proposition 13.2(1 

\Im,d( z )\ < 27rC m , d e- /312 for any < ft < 1 - p, 

< 27r(7 m>d e- /32Z for any < ft < P- (5.5) 



Term Z ajb — (if) a , P u i^ b ). Using the expression for Z ajb , see (13.24(1 . we obtain 
{^{t-K)^ b )-(^PM 

" a + 6vvl + 2al + 267 vl-4a^ 6 v 1 ' 

(|5.6() is analytic for a, 6 G (—1/2, 1/2) because the two terms in the bracket 
are 1 + 0(a + b) when b + a — > 0. 

Term (ip a , P u K m ^ipb) ■ Using (|3.15() and (|3.17() one obtains 

(i> a ,PuK m4 i> b )=Z(-b) f°d* f dy ImAx ~. y) e~^ + ^ e y^- b \ 

Ju J-oo 

(5.7) 

where the function I m ,d(z) is given in (|3.18j) . Z(—b) is analytic if b > —1/2. 
Thus the integrand is bounded by C mjd e~ x ^ 1+a+p ~2 S >e y ( /3l ~ b+p ~2\ The condi- 
tion f3i < 1 — p implies that the integrand is exponentially decaying in x — y 
provided that a > — 1/2 and b < 1/2. Thus r.h.s. of ()5.7j) is real analytic for 
a, b G (-1/2,1/2). 

Term (i/; a , Q u (t - K m)d )P u (t - P u K mid P u ) _1 P u (l - K m>d )i/; b ) . The object to 
consider is 

poo poo 

dx dyf a (x)(l-P u K m4 P u )- l (x,y)g b (y) (5.8) 

J u J u 

with f a (x) = -(K* md Q u i) a )(x) and g b (y) = ((1 - K mid )i/; b )(y). Explicitly, 
fj x ) = e -<i-P>» f°d™/"° ^-^ ^M^-^lUg^) 



—oo 




-^(a)-^-^-^ y d^e-^a+P-l) 7 ^.^ z ) (5.9) 
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and 

9b ( x ) = e~ bx - e-^ x Z(-b) f dz Im4 ^~ z) e z ^- b ). (5.10) 

J -oo —1t\% 

Using (1 - P u K mA P u y x = 1 + (1 - P u K mA P u y x P u K mA P u we rewrite Q 
as 

(fa,Pu9b) + (fa,P u (l ~ P u K m4 P u y l P u ~g b ) (5.11) 
with g b = K mA P u g b . Using ()5.5|1 one deduces that, for a G (—1/2, 1/2), 

/„ G L 2 ((0, oo), e^dx) for all p < 1/2 (5.12) 

and, for b G (-1/2,1/2), 

g b G L 2 ((0, oo), e-^dx) for all p < -b < -. (5.13) 
This implies that 

POO 

\(f a , P u 9b)\ < C dxe-^ 2+1 ^ p)x (e- bx + e- {l3l+p -^ )x ) (5.14) 



for some finite constant C . One has P2 + \ — p + b — >| + 6>0as/9 2 -^p 
and /3i + /?2 > 0, from which it follows that (f a , P u gb) is real analytic for 
a, 6 G (-1/2,1/2). 

For the other term in 1)5.11)1 we have to compute g b (x). We obtain 

# fe (x) = e (p ^ dw dz '2 2jr . 2jri -e w(2 p) g b (w). 

(5.15) 

Then 

POO pO 

\~g b (x)\<Ce- {Pl+p --i )x dw dze^ 1+f32)z (e- {b+ ^-p) w + e -^ + ^) 

J u J —00 

(5.16) 

for some constant C. One has (3\+ p — | — > | as /3j — > 1 — p, /3 2 + | — p — 6^ 
I — 6 > as /3 2 — ► P, and /3i + /3 2 . It follows 

g b G L 2 ((0, 00), e px dx) for all p < 1/2. (5.17) 

From this one deduces that (/„, P u (l — P u K rnA P u )~ l P u g b ) is also real analytic 
for a,b E (—1/2, 1/2). 
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A second representation. Another way of eliminating the divergence consists 
in subtracting the quantity (ip a , (1 — K md )P u ip b ) from each term. This leads 
to 

r.h.s. of Q = (Z a>b -(ilj a ,P u il> b )) + {iP a ,K mjd P u iP b ) (5.18) 

-<Va, (1 - K m4 )P u {t - PuK^dP^P^l - K rn4 )Q u ^ b ). 

By the same argument as above one shows that the terms in this second 
representations are analytic for a,b G (—1/2, 1/2). □ 

By Propositons l5.1l and l5.2l one can take the limit b — > —a without loosing 
the identity (|4.1Uj) . We remark that lim^_ a (a + b)Z a ^ b = 1. In addition, by 

2(2G? — 777, 

lim {ip a , (1 - K m4 )^ b ) - P u ^ b ) = u + -. (5.19) 

b^-a i — a Z 

Let us denote the limit 

lim G a '\u) =G {u), (5.20) 

b^— a 

with the a-dependence understood implicitly. Then by ()5.4j) 

2ad — m 



/ Ann 777 \ 

G (u) = lim G a > b (u)= (u + - r )+{t/>a,PuKmM>-*) ( 5 - 21 ) 

b^-a V ± — or / 

-(Va, Qu{t ~ K m4 )P u {l - PuK^P^P^l - K m4 )^_ a ). 
Alternatively, using ()5.18|) one can write 

/ 2/7/7 777 \ 

G Q (u) = lim G a \u)= (u + - r )+(iP a ,K mid P u ^ a ) (5.22) 

b^-a V ~ — a 1 / 

-(i/) a , (1 - K m4 )P u {l - P u K mA P u )- l P u (l - K m4 )Q u ij_ a ). 
With the same convention, let us denote 

lim Pq' 6 = P . (5.23) 

b^—a 

Recall that by construction Po is the probability measure for the family 
w(i,j), i,j G N, of independent exponentially distributed random variables 
such that w(0, 0) = 0, w(i,0) has mean (1 — p) _1 for i > 1, w(0,j) has mean 
for j > 1, and otherwise w(i,j) has mean 1. This is precisely the w- 
marginal of (Q in Section El Also, by definition, h°(j,r) = ho(j,r) pathwise. 
Thus one concludes 
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Proposition 5.3. For < C\ < c 2 < oo it holds 

[ 2 duQ({h°( 3l T)<u}) = f 2 duP ({h ( 3l r)<u}) 

J CI J CI 

= F(c 2 )G (c 2 )-F(c 1 )G (ci). (5.24) 

To prove Theorem l2.4l one has to investigate the asymptotics of F{u)Gq(u) 
under the scaling (j2.12|) . 

6 Edge scaling 

Following ()2.12|) we set 

2m = t - 1 = [(1 - 2x)t + 2w(l - 2p)x 1/3 t 2/3 - sx 2/3 t 1/3 \ , 

2d = j = L(l - 2p)t + 2^ X 1/3 t 2/3 J , (6.1) 
u = [t + sx-^l 

with x — p(l — P)- The, by Proposition 15. ^ the proof of the limit (j2.14j) in 
Theorem 12.41 reduces to the large t limit of 

X-'/H^Goit + sx-^t 1 / 3 ) (6.2) 

where the prefactor takes into account the scaling of as well as the large 
t limit of the Fredholm determinant 

det(l - P u K m4 P u ) (6.3) 

on L 2 (R + ). This latter limit has been studied by Johansson, see Theorem 

i.6 in ma. 

Theorem 6.1. (Johansson) Let u = t + sx^ 1 ^ 3 ^^ 3 and m,d as in \6.1\) . 

Then 

lim det(l - P u K m4 P u ) = F GVE (s + w 2 ), (6.4) 

t— >oo 

where Fque is the GUE Tracy- Widom distribution function \2J$ . 

In order to state the limit of Go we have to introduce some auxiliary 
quantities. We define the functions 

<p w , a (z) = Ai(w 2 + s + z ) e ^ 2 +^) e -^ 3 ( 6 .5) 
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and 



S w , s = s+ dxdytp WjS (x + y), 

= e< / dycp_ WjS {y + C)(l- I dx ( p W;S {x + y)\ (6.6) 

</ 1R-|_ J ]R+ 

= e-^fl- / dj^ w>a (a; + 0). 
These functions can be written as a single integral using the identity |D] 



1- / dx</? w , s (x + y) = I dxip WjS (x + y) (6.7) 

for io > (for w = the same holds but only as a improper Riemann 
integral). Using a contour integral representation, in case w > 0, one rewrites 
S 1 ,,, o as (see Section ID. 3D 



S WjS = / dxdytp w>s (x + y). (6.8) 
Jr 2 _ 

It is easy to see, using the super-exponential decay of the Airy function, that 
® w ,s € L 2 (R + ) for w > 0, \f^ )S e £ 2 (R+) for w > 0, and in the case w = 0, 
*&o,s - 1 G L 2 (R+). Finally we denote by K^ q the operator with kernel 

AAi,,(£i,6)= / dzAi(g + z + £i)Ai(g + * + (6.9) 

*/R+ 

i^Ai,g is the Airy kernel shifted by g. 
Theorem 6.2. 

lim x 1/3 t~ 1/3 G (t + sx" 1/3 t 1/3 ) (6.10) 



t— >oo 

= <S'|-!il|,S 4" 



2 



= #(s + w ,w). 

Theorem 16. II and Theorem 16 .2| in conjunction with Proposition 15.31 fur- 
nish the proof of Theorem 12.41 and hence of our main result Theorem 11.21 
Beyond the existence of the limit, they also implies Corollary 11.31 

Proof of Theorem \6.2L 

A change of variable: from t to N. We change variables with the effect 
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to have N instead of t as large parameter. This is not really necessary but 
simplifies our computations. Let us define N = m — d = (r — 1 — j)/2 and 
N + a = m + d. The relevant parameters for what follows are a = 2d = j 
and u which, in terms of N, are given by 



p 2 p 2 



It 



+ (|^ 2 (l - P) + »(1 - 2p)) (1 ^ /3 iV'/3 + (i), (e.ii) 

4 + 2m < 1 -fV° + (V + sf-fV* + 0( i). 

p z p^ V3 / p z 

Moreover the scaling x^ 1 / 3 ^/ 3 writes as 

x -l/3 f l/3 = rJV 1/3 + (Q_ 12 ) 

with 

K = p- 1 (l-p)-V3_ ( 613) 

After edge scaling the terms of Go will be expressed via the functions Hn 
and if at, defined as 

k,N 1/3 

H N {y) = Zi^—^—I^^iu + ynN 1 / 3 ), (6.14) 
nN 1/3 

H N {y) = Z(a)- l ^—I N+aM2 (u + yKN 1/3 ). 

Using the bounds on / and I of Section [7| we obtain, for any (5 > fixed, 

\H N (y)\ < Cpe-^, (6.15) 
\H N (y)\ < Cpe-^y 

for some Cp > independent of N and y > 0. Moreover we also have the 
point wise convergence 

lim ffjvfo) = Ai^ + s + ^e^+^e-^ 3 = <p w M, (6-16) 
lim tf^y) = Ai(w 2 + s + y)e- w ^ 2+s+ ^e 1 s w3 =ip^ s (y). 

N^oo 

We simplify the notations in this proof by setting K = iGv+a/2,a/2- 
The estimate of the terms for «;~ 1 iV -1 / 3 G ! o(it + skN 1 / 3 ). 
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First term of \5.21)) and \5.22}) . Using (|6.11J1 . 2d = a, and m — d = N, we 

have 

2ad — m 

lim u H — s — = s. (6.17) 



N^oo J — a 



Second term of \5.21\) for w > . We compute the limit 

lim K^N^i^PuK^a). (6.18) 

N—*oo 

ip- a is eigenfunction of R, see (13. 15 J) . Thus we have 

poo poo 

(ij> a , P u KiP_ a ) = Z(a){iP a , P tt LP_V_.) = rf/ 3 / dx / d^(x + y). 

Jo Jo 

(6.19) 

Using the bound (|6.15|) we can apply dominated convergence. Then by the 
pointwise limit ()6.16|) we have 

poo poo 

lim K- 1 N- 1 / 3 {i> a ,P u K^ a )= / dx dy<p w , s {x + y). (6.20) 
^- +oc Jo Jo 

Second term of \5.22}) for w < 0. This case is analogous to the previous one. 
One obtains 

poo poo 

]im K- 1 N- 1 ' 3 U> a ,KP u il>- a )= dx dy<p- WtS (x + y). (6.21) 

Jo Jo 

The sum of fl6.17|l and (j6.2()jl . resp. (|6.21jl . yields S\ w \ tS as the first term in 
Theorem 16.21 

Third term of \5.21\) for w > 0. The third term of n~ 1 N~ 1 / 3 Go(u), including 
the pref actor —1, is 

k^N-W^n^nVn) (6.22) 

with 

$ N (x) = (K*Q u i> a )(x), 

*N(y) = ((l-K)i>- a )(y), (6.23) 
A N (x,y) = {P^l-PuKP^P^x.y). 

To establish the scaling limit, x = u + ^kN 1 ^ 3 , we define the rescaled quan- 
tities 

i>L a (0 = ^-aiu + ^N 1 ^), (6-24) 
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where 

= e -a(u+^N^) e - W ^ (6>25) 

Then fid 22 j) becomes 

($ N ,A N V N ) (6.26) 

with 

= (^UV4)(6), 

<M6) = (6.27) 

We can rewrite the <3>at, ^at, and Kj^ by using the functions H N and H N , 
<Md) = e^f/ dyH N (y + ^)- [ dxdyH N (x + y)H N (y + frj) , 
<M&) = e~ w ^(l- [ dyH N (y + &)), (6.28) 
= e-^e^ 2 / dz# A r(:r + £ 1 )# JV (x + £ 2 ). 

./R+ 

We want to avoid to write always the projection P Q . Therefore from now 
on (•, •) refers to the scalar product in L 2 (R + , dx), \\ ■ || is the corresponding 
norm, and the integral operators act in L 2 (H + , dx). First let us consider 
w > 0. Let us denote A = (1 — K^ w 2 +s y 1 . Then for finite N we have the 
bound 

\($ N ,A N * N ) - ($ w>s ,Ay WjS )\ < \\$ N \\ \\A N -A\\ \\V N \\ (6.29) 
+\\$ N - QuJ \\A\\ \\V N \\ + \\® W J \\A\\ \\V N - V WtS \\. 

In Lemma IFTTTT1 we will prove that $at converges to $ WtS , converges to 
ty w , s , and An converges to (1 — K\ itW 2 +s )~ 1 in operator norm (in L 2 (R + ) 
according to our convention). This implies 

lim ($ N , A N V N ) = (<S> W>S , (1 - K AijW 2 +s )- 1 ^f WjS ) (6.30) 

N— >oo 

which is precisely the last term in (jfi.lOj) . 

For the case w = we have to modify slightly the argument. In this case 
^ N is not in L 2 (R + ), but - 1 = - L dyH N (y + L 2 (R + ) and 

converges to — J R dx(po tS (x + £) in the N — > oo limit. We write 

($iv, A^v^at) = ($ N , A N (V N - 1)) + ($jv, Ajvl). (6.31) 
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Here 1 denotes the constant function l(x) = 1, for all x 6 R. By the same 
argument as for w > 0, the first term in (|6.31|) converges as 

lim ($ N , A N {* N - 1)) = ($ w , s , (1 - Km^+s)- 1 ^ - 1)). (6.32) 
The second term in (I6.31j) can be rewritten as 

($^,^1) = / d£id£2$tf(6)4v(ei,6). (6.33) 

In Lemma [6.41 we will prove the convergence 

lim ($ N , A N 1) = ($ WjS , (1 - K Aiiw2+ ,)- 1 l). (6.34) 

Third term of \5.22]) for w < 0. The computations for this case are as before 
and the third term of ()5.22|) converges to 

(*_ w , s , (1 - K MtV fl +a )- 1 ^- W)t ). (6.35) 

Since K^ w 2 +S is symmetric, this concludes the proof of Theorem 16.21 

□ 

Lemma 6.3. Let w > 0. Then 

lim H^jv - $ MjS || = 0, 

lim H^-^ll =0, (6.36) 

N— >oo 

lim \\A N - (1 - -^Ai,«;2 +S ) _1 || = 0, 

N— >OC 



where the functions $at, ^jv; ^ e integral kernel A^ are defined in \6.21 ) 
and |Og) . 

Proof. Convergence of <&n- Let us consider, for any fixed £ > 0, the function 
$jv(£) defined in (JH2EJ). We first show that 

lim = $^ s (£) (6.37) 

pointwise. Using the exponential decay ()6.15jl of Hn and if at, we apply 
dominated convergence and exchange the integrals with the iV — > oo limit. 
Then using the pointwise limit of H^, see (16.16)1 . one obtains 



lim = e™* / dy<p_ W)S (y + £)-e w t dydx<p_ w>s (y + £)(p w , a {x + y), 

(6.38) 
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which is precisely $ WtS (£)- By the exponential decay (j6.15|) with ft > w 
it follows that $jv 6 L 2 (R + ). Moreover, & N is uniformly bounded by an 
integrable function since, for all ft > w, 

IMOI < Cje'^e-^e-^-^. (6.39) 

Therefore, by dominated convergence and by pointwise convergence (jfi.lSSjl . 
we obtain 



lim ||$tv - $w,s\\ 2 



[ dC lim |^(0-^,.(0r = 0- (6- 4 °) 
Jr+ n ~*°° 



Convergence of^N- Let us consider, for £ > 0, the function ^n(0 defined 
in JS2H|. We first show that 

lim = * w ,-(0 (6-41) 

pointwise. As before, f)6.15j) allows us to exchange the limit and the integral 
with the result 

lim tfv(0 = (l - /°° dj/^y + 0) • (6.42) 

iV^oo ^ Jo J 

Then by (JHHEJ), with > w, G £ 2 (1R+) and ^ is bounded by an 
integrable function. Therefore 

lim ||^v - V w ,s\\ 2 = [ d£ lim |^v(0 - ^(OT = 0- (6.43) 



'R+ 

,2 



Convergence of Am- Denote q — w + s. We want to show that 

lim ||(1 - K r N y l - (1 - Kax^W = (6.44) 

N— >oo 

in operator norm. Assume that we can show: 

1) liniA^oo \\K r N - KmA = 0, 

2) \\{t - Ka^)- 1 ]] < oo. 

Then with the notation K £ = K^ q — K r N we have 

w(i - K^r 1 - (i - ka^w 

= ||[(1 + (1 - K M , q )- x K E )- 1 - 1](1 - Ka^W 

< 11(1 - ^Ai,,)" 1 !! ^ ||(1 - K^rTWKeir (6.45) 

n>l 
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which converges to as N — ► oo. So we have to establish properties 1) and 
2). Let us start with 1). For fixed x, y > 0, using (j6.15|) with (3 > w, K r N (x, y) 
is uniformly bounded by a function which is integrable in R+. Therefore, by 
dominated convergence and the pointwise limit 

lim H N (x + z)H N (y + z) = Ai(q + x + z) Ai(q + y + z), (6.46) 
one obtains 

lim \K r N (x,y)-K M>q (x,y)\ (6.47) 

/■oo 

< / dz[ lim H N (x + z)H N (y + z) - ki(q + x + z) ki(q + y + z)} = 0. 
Jo N ^°° 

Moreover, it is easy to see that K r N is a Hilbert-Schmidt operator with norm 
uniformly bounded in N, and so is K Aiq . Therefore 

lim \\K r N (x,y) - K Mq (x,y)\\ 2 < lim \\K r N (x, y) - K M>q (x, y) ||| s 
= [ dxdy lim \K r N {x,y) - K M>q {x, y)\ 2 = 0. (6.48) 

Next consider point 2). 1 — K^i^q is invertible as bounded operator, since for 
every g 6 R, ||-^Ai,g|| < 1- To establish the claim, let us denote by K Ai the 
standard Airy operator with Airy kernel K^q. Then 

\\K Aijq \\ L 2 (u+) = \\K M \\i?{[q,oo)y (6-49) 

is an operator on L 2 ([q, oo)) which is Hilbert-Schmidt. Therefore the 
norm of -f^Ai,? on L 2 ([q, oo)) equals its largest eigenvalue, Xo(q). In [24 j it is 
shown that \o{<l) is monotonically decreasing in q and is strictly less than 1 
for any q > — oo (it converges to 1 as q — ► — oo). □ 

Lemma 6.4. Let w = 0, then 

lim ($ N , A N 1) = ($„,., (1 - ^Ai,.) -1 !)- (6.50) 
Proof. We first rewrite 

(* N , A N 1) = ($ N , 1) + ($ N , (1 - K r N )~ l K r N l) . (6.51) 



The first term, ($^,1) = f dx^jy(x), converges to L dx$ ,s 



. x 



($o,s, 1) as iV — > oo because $n G L 1 (]R + ) with norm L 1 (R + ) uniformly 



bounded in N. For the second term, define $/v = (1 — K^*) 1 $\r. From 
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Lemma f6. 31 it follows that, for N large enough, <&/v e L 2 (R + ) with uniformly 
bounded norm. Then 

($jv, (1 - K^K^l) = [ dx [ dy$ N (y)K r N (y,x) = [ dxD N (x). 

(6.52) 

Dn(x) is bounded by an L 1 function independent of N. In fact, using the 
representation ()6.28|) of K r Nl one writes 



\D N (x)\ < / dX\H N (x + X)\ dy\H N (y + \)\\$ N (y)\ (6.53) 

JR+ JTR.+ 

< [ dX\H N (x + X)\( [ dy\H N (y + X)\ 2 ) 1/2 \\$ N \\ 2 . 
Jr+ y Jn+ ' 

Since Hn and decay exponentially, see (j6.15j) . we have \D^(x)\ < Ce~ x 
for some C > independent of N. Therefore by dominated convergence we 
obtain 

\im(<S> N ,(l-K r N )- 1 K r N l)= [ di Jim / dy$ N (y)K* N (y, x). (6.54) 



This last integral can be interpreted as the L 2 (R + ) scalar product between 
$at and K r N {-,x). By Lemma l?TTTT| $jv converges to (1 — Kai, s )~ 1 ^o,s- More- 
over, K r N (-,x) converges to Kai,s(-,x), thus 

lim (Q Ni (1 - K r N )- l K^ N l) 

= [ dx [ dy(l-K Ai ^-%, s (y)# Ai)S (y,x) (6.55) 
Adding the first and second term one obtains (|6.5()jl . □ 

7 Asymptotics used in Section M 

In this section a and u are defined as in (|6.11j) . First we summarize the 
asymptotic results required. Let m = N + a/2 and d = a/2. 

7.1 Asymptotics of I m ,d{u + yhiN 1/3 ), < p < 1 
Oi l For fixed j/6R, 

Z(a)nN 1/3 I md (u+ynN 1/3 ) = -2m Ai(w 2 +s+y)e w(w2+s+y) e~^ w3 +0(N- 1/3 ). 

(7.1) 
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17.11 2 For < L < ?/ < eN 2 / 3 , e > small enough, L > large enough, 

\Z(a)KN l / 3 I m4 (u + |/KiV 1/3 )| < Ce"^ 3/2 (7.2) 

for some C > 0. 

17711 3 For y > eiV 2 / 3 , e as in ()7.2|) . 

^(a)/^ 1 / 3 /^ + ynN l '*)\ < Ce~^ £l/2yNl/3 (7.3) 

for some C > 0. 

7.2 Asymptotics of 7 m ,d(w + ijkN 1/s ), < p < 1 
17721 1 For fixed t/6R, 

Z{a)KN 1 ^ 3 L m d^u + ynN 1 ^ 3 ) = 2ni Ai(w 2 + s + y)e- w{w2+s+y) e^ w3 + 0(N- 1/3 ). 

(7.4) 

I7.2L 2 For < L < ?/ < eN 2 ^ 3 , e > small enough, L > large enough, 

|Z(a)«iVV 3 / m>d ( u + yKN l ' 3 )\ < Ce~-y /2 (7.5) 

for some C > 0. 

17721 3 For y > eiV 2 / 3 , e as in (I73|) . 

^(a)/^ 1 / 3 /^ + i/^iV 1 / 3 )! < Ce-^ 1/2j/iVl/3 (7.6) 

for some C > 0. 

Proof of \7.1\j . We have to estimate 

^ri_ p 

with 

f N ( z ) = - Z {u/N + ynN~ 2 / 3 ) + \n(z + p) - (1 + a/iV) ln(l - p - z) (7.8) 

for any fixed t/6R. 
Let us define 

f oc (z)= lim /j^z) = + + p) ~ (1 ~/ )2 ln (! ~ P ~ *) ( 7 - 9 ) 

31 



Figure 3: Integration path used for the asymptotics for fixed y. 



and find a steep descent path for it which is close to the steepest descent one 
for z close to the critical point, which is the solution of 

d/ooO) 1 1 (1 - pf 1 , 
az fr z + p p z 1 — p — z 

There is a double solution for z = 0, thus f^(z = 0) = 0. Moreover, 

d 3 /oo(2 



dz 3 



=o p 3 (l-p)' 



(7.11) 



Therefore we choose as integration path the one shown in Figure El The 
chosen path is a steep descent path for /^ as is discussed now. 

The path j 2 is given by {z = te~ l7T ^ 3 ,t G [0,2(1 — p)]}. The real part of 
/do on 72 is then 

Re(U = ~ + J In ((„ + i*) 2 + J* 2 ) " In (d - P - \tf + \f) . 

(7.12) 

Therefore 

di^e(/ 00 ) t 2 (2p(l-p)+t(l-2p)+t 2 ) 

dt 2p 2 (p 2 + pt + t 2 )((l-p-t) 2 + t(l-p))' l ' ] 

The denominator is positive and it is easy to see that the numerator is always 
strictly positive for t G (0, 2(1 — p)] and for all p G (0, 1). Therefore 72 is a 
steep descent path, and by symmetry 71 is a steep ascent path. 
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The path 73 is given by {z = 1 - p + V3(l - p)e i</3 , <p G [— 7r/2, tt/2]}. On 73 

Jfe(/oo) = - (1 " C ° S ^ + 5 ln(l + 3(1 - pf + 2^3(1 - p) cos <p) 

p z 2 

-ii^l!ln(V3(l-p)), (7.14) 
P 

which implies 

diteC/oo) (l-p)V3 , VS(1- P ) 



dcos<p p 2 1 + 3(1 - p) 2 + 2^(1 - p) COS if 

(1 - p)v/3(2y/3(l - p) cosy? + 3(1 - p) 2 + 1 - p 2 ) 
p 2 (l + 3(1 - p) 2 + 2>/3(l - p) cosy?) 



< 0. 



Thus the path of Figure El is a steep descent path 

The first consequence is the following. Denote by r 1 _ p (S) = Ti_ p \,. <s 
the part of the path Ti_ p closer than S to the origin. Then for any 5 > and 
N large enough, 



< e NfN{0) O{e-^ N ) (7.15) 



for some p = p(£) ~ 5 3 > 0. Remark that 

e iv/ w (o) = Z ( a y\ ( 7 16 ) 

Consequently we need to estimate the integral close to z = on Ti_ p (5) only. 
We use the Taylor expansion, 

f N (z) = fN(0) + f N (0)z + ^(0)z 2 + lf^(0)z 3 (7.17) 

+0{\z\* max \fP(z)\). (7.18) 

0<|z|<o 

Some computations yield 

f' N (0) = -N- 2 '\y + s)K + 0{N- 1 ), 

/£(0) = N~^2wk 2 + 0(N~ 2 / 3 ), (7.19) 

/#(0) = 2k 3 + 0(N- 1 / 3 ), 



and \f^ v (z)\ = 0(1) for \z\ < 5. The change of variable r = N 1 I 3 kz leads 
to 

jV/y(z) = Nf N (0) - r(y + s) + wr 2 + ^r 3 + 0(r, r 4 )^- 1 / 3 . (7.20) 
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Consequently 

e -Nf N (0) f e NMz) dz = K -l N -l/3 f dre -r(y+s) +W r^ 

+ K^N- 1 / 3 [ dre ^+-)W+|T3/0(r^)JV-V3_ 1 Y (7 _ 21) 

The last term can be estimated using that \e x — 1| < l^e'^', i.e., using 

The term in the exponent is of the form — r(y + s)xi + wt 2 X2 + | r3 X3 f° r 
some Xi, X2, Xz- By taking <5 small enough xii X2, X3 can be made as close to 
1 as desired. Thus the second integral converges and the error term in (j7.21j) 
is of order 0(N~ 2 / 3 ). 

Finally we estimate the leading term 

K-ljV- 1 / 3 / dre -r( y+S )W + ir3_ ( 72 ^ 

J KNV 3 ri- p {5) 

Deforming the integration path from r(5)A^ 1 / 3 K to T(5)N 1 ^ 3 k— w, one obtains 

e (y+s)w e l^ K -l N -i/3 I dTe - T{y+s+w 2 )+ i T 3 ( ^ 24) 

Jf(s)n 1 / s k 

up to an 0(e~^ N ) error. By extending the integral to e ±l7T ^ 3 oo one picks up 
an error at most 0{e~^ lN ) ) again. But 

/ dze^~ xz = -2iriAi{x), (7.25) 

where T(oo) is the path joining with e ±l7r ' 3 oo by straight lines oriented with 
imaginary part decreasing. Note that in ()D.2|) the orientation is the opposite. 
The error term 0(e~^ N ) can be bounded by 0(N^ 1 ^ 3 ). Thus putting all the 
terms together we have proved that for any fixed y 

lN +a /2,a/2(u + ynN 1/3 ) Z (a) kN 1/3 = -2m Ai(w 2 + s + y) e ^+^+y) e ~W 

+ 0{N~ 1/3 ). (7.26) 

□ 

Proof of [71$ . Let us define y = yuN^ 3 G [LkN~^ 3 , we], 

f N (z) = -z(u/N + y) + ln(p + z) - (1 + a/N) ln(l -p-z), (7.27) 
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7i 



72 



Figure 4: Integration path used for the asymptotics for large value of y. 72 
is the dashed line. 



and 

g(z) = -z(p- 2 + y)+ ln(p + z)- (1 ~ P) ln (! ~P~z)- (^.28) 

P 2 

g has a real positive critical point at 

z^y^n-^ + Oiy). (7.29) 

Let = y 1 ^ 2 ^ 3 ^ 2 . Then as integration path we choose Ti_ p = 71 A72, where 
71 = {^|z = Zq — it,t E [v3z , \/3^o]}) 72 the path used in the case of fixed 
y restricted to Re(z) > z , see Figure |H Since y > 0, the path 72 is steep 
descent because it is for y = and Re{z) > on 72. Thus we only have to 
check it on 71. On 71 we have 

Re(f N (z)) = -, (^+^+Iln((p+ 2o ) 2 +t 2 )-i(l + ^) ln((l-p-^ ) 2 +t 2 ) 

(7.30) 

and 

dRe(f N ) _ (t 2 + (p + z ) 2 )a/N + 2p-l + 2z 

dt ((p + z )i + ti)((l-p-z o y + ti)- 1 • ] 

The denominator is obviously positive. 
Next consider the numerator 

M = (t 2 + (p + z ) 2 )a/N + 2p-l + 2z . (7.32) 

For p G (0, 1/2), a/N = (1 - 2p)/p 2 + 0(N^ 1 ^ 3 ), M can be rewritten as 

M = a/Nt 2 + (1 - 2/>) ((1 + Vp) 2 - 1) + 2z + (^(iV™ 1 / 3 ) > 
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for N large enough. For p E (1/2, 1), a = (1 - 2p)p 2 + 0{wN~ 1 ' 3 ) and M 
can be rewritten as 



M = ^-t 2 + 2z Q 1 —^ + ^ 2 ^-^ + 0{wN-^) 
JS p p z 



Since z = y 1/2 K~ 3/2 > L 1 / 2 k~ 1 N~ 1 /' 3 , and that z < e^ 1 < 1, both the 
quadratic term in z and the 0(wN^ 1 ^ 3 ) are dominated by 2z (l — p)j p for 
L > 1, e < 1, and iV large enough. Thus M(t = 0) > L 1 / 2 /^ 1 ^ 1 / 3 > 
for L large enough. By monotonicity of M in t we have to check that M(t = 
V^o) > 0, the maximal value which t takes in 71. But M{t = V3zo) — M(t = 
0) ~ Zq, which is dominated by the linear term in zq again. Thus M(t) > 
for N large enough and for all t E [— V3z , V3z ]. 

We have shown that for all p E (0, 1), dRe ^ N) < for t > 0, dRe ^ N) > for 
t < for L 1, e -C 1, and A/" large enough. Thus 71 is a steep descent path. 
Therefore, if we denote by Ti^ p (5) c the portion of with \z — zq\ > 5, 



dze NMz) 
ri- P (5) c 



for some p > 0. 

Finally we have to evaluate the contribution coming from Ti_ p (S), the 
portion of ^i- p with |z — z \ < 5. The contribution of the part in ji is 
estimated as follows. On ji, 

Re(f N (z)) = f N (z ) - jftizo) + 0(t 4 ). (7.34) 

Some computations leads to /^(^ ) = 0(wN~ 1/3 ) + (2k 3/2 + 
0(wN~ 1/3 ))y 1/2 + 0(y). For L > 1, e < 1, and iV large enough, it fol- 
lows that f'^(zo) > I/? 3 / 2 ?/ 1 / 2 . On the other hand, the term proportional 
to t 4 is much smaller than the t 2 term. In fact, for < t < V^Zq ~ 2/ 1 / 2 , 
i 2 < yC(l) < y l l 2 0{e l l 2 ) < /^(^o). Therefore 



dze NfN{z) 



71 



= e 7V/,v(20) ?/" 1/4 ^~ 1/3 O(l). (7.35) 



In fl7.33|l and ()7.35j) we still have to evaluate /at(^o) — /at(0). A compu- 
tation leads to 

A7iv(*o) - Nf N (0) = -sy l ' 2 xi + wy X 2 - (7.36) 



36 



for some Xi, Xii X3 which can be made as close to 1 as wanted by choosing e 
small enough. On the other hand, since we have y > L, for L large enough, 

Nf N (z )-Nf N (0)<-±y 3 /\ (7.37) 

Thus the contributions ()7.33|) and ()7.35|) can be bounded by 
TO < Z(a)- 1 e-^e-^ 3/2 C(l), 

(TT351) < Z(a)- 1 iV- 1 / 3 e-^ 3/2 L- 1 / 4 C(l). (7.38) 

The final step is to bound the contribution coming form 72 U Ti- P (5). 
From the asymptotics of fixed y one has, using z = te t7T ^ 3 , 

Re(f N (z)) = f^-ty + sKN-^ + OitN- 1 ) (7.39) 

+W ^N- 1 ^- + 0{t 2 N- 2 ' z ) - -t 3 + 0{eN- x '\t A ). 
2 3 

In this case, the parameter t takes values in < 2z < t < 5/y/3 <C 1. 
Moreover recall that z > L 1 / 2 k~ 1 N~ 1 / 3 . In the term linear in t, y dominates 
the others for large L. For the minimal value taken by t, the quadratic term 
is ~ LN' 1 , and the cubic term is ~ L^^N" 1 . Thus for large L, the cubic 
term dominates the quadratic one. But since t < 5/v3, the quartic term is 
also dominated by the cubic one. Therefore, 

Re{f N {z)) < MO) - — - — < MO) - ^ - — . (7.40) 

Thus 

/ d^e^W < 2e ^(°)e-^ 3/2 / dte~^ 3iV < ^ f ' O(l). 
A 2 ur 1 _,( ( 5) Jo ^ 1/3 

(7.41) 

From ()7.38|) and (|7.41|) the desired bound follows. 

□ 

Proof of \7.3j) . The proof of this bound is based on the estimate ()7.2|) . We 
use (17. 2 j) for y = ek/2 with the result 

M*) = < /*(z)| g ^ /2 -2- 3 / 2 v%iV- 2/3 , (7.42) 



because 2 > \Js/2k 1 and y — en/2 > y/2. It follows that 



W/2,a/ 2 (« + yKN l / z )KN l ' z Z{a) 



< Ce~ £Z/2N,2 e~^ yNl/A,A < Ce~^ yNl/3/4 . 

(7.43) 
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□ 

Proof of H7.4\ ), and J7.6| ). The asymptotics of I are similar to the 

one of /. Instead of computing everything again, we show that, via the 
transformation p i— > 1 — p and w i— > —w, one obtains essentially the same 
integrals as already studied for /. More precisely, we have to estimate the 
asymptotics of 

W(N,p,w) = Z{a)KN x l*I N+aM2 (u + yKN x l*) (7.44) 



N 

d^e" 2(u(7V ' p ' w)+s/K(p)iVl/3) K(p)A^ 1/3 



p 



1 - T , 

1-/0 y 



N+a(N,p,w) 



and 



W(N,p,w) = Z{a)- l KN l l*I N+a/%a/2 {u + yKN x l*) (7.45) 

/ \ N+a(N,p,w) 



N 

Here the dependence in N,p,w of u,a,K is displayed explicitly, since it is 
needed below. A simple calculation shows that 

W{N, 1 - p, -w) (7.46) 

/ \ N+a(N,l-p,-w) 

(l- £ ) 

d 2e -*(«(^.i-p-w)-HfK(i-p)iv 1/3 ) K (i _ p)ArV3 V 'V 



iV 



Let us define M = N + a(N, 1 — p, —w), then with an explicit but lengthy 
computations one establishes 

u(N,l-p,-w) = u(M,p,w) + 0(l), 

K{l-p)N 1/3 = /<p)M 1/3 + 0(l), (7.47) 
N = M + a(M, p, w) + 0(1). 

(|7.47jl implies that in the asymptotics of W(N, 1 — p, —w) are the same as 
the asymptotics of —W(M,p,w), since the corrections of 0(1) in (|7.47jl are 
negligible for the asymptotics ()7.1|) . ()7.2j) . and (|7.3|) . 

□ 
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A Scaling functions 



The scaling function g(s,w) denned in (jl.lSj) is precisely the one derived in 
(jfi.lOJ) except for the shift to s i— ► s — w 2 , since in (Ifi.K)j) g(s + w 2 ,w) is 
obtained. Below we establish that g is identical to the Baik-Rains scaling 
function g BR [2] g(s,w) is continuous at w = and even in w. <?br(s, w) 
has the same properties. Thus it suffices to consider w > 0. We first rewrite 
(jfi.lOJ) in another form by moving the factor s from the integrand to the limit 
of the integral, since later on we have to deal with the derivatives -g^g(s,w). 
Define the functions 



$ a (x) = / dyAi(x + y) / dze wz M(y + z), 

JO J -oo 

= f dye wy Ai(x + y), 
p s (x,y) = (l-P s K Ai P s )-\x,y). 
Then by using the representation ()6.8|) for S W;S , ()6.10p rewrites as 

f° 

dx / dy Ai(x + y)e w(x+v) 

I J —oo 

POO 

+ I dx dy§ s {x)p s {x,y)^ w {y) 



(A.l) 



g(s,w) 



-±w 3 

e 3 



The relations with the functions of (|6.10|) are $ s 



x 



[x 



u2 (x - s), and y) = (1 - P -^Ai, s -Po 



u2 (x - 
- s,y 



(A.2) 

a), 
s). 



Baik and Rains 0, see also [IB], use the Riemann-Hilbert techniques 
and arrive at a limit function, <7br> which is given as the solution of a set of 
differential equations. More precisely, with a = a(s, w), b = b(s, w), q = q(s), 
one considers 



—a = qb, 

OS 

—b = qa - 

as 



wb, 



(A.3) 



and 



d 

—a = q 2 a - (q + wq)b, 
d 

7—6 = (q - wq)a + (w 2 
Ow 



q 2 )b. 



(A.4) 
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Here q = q(s) is the Hastings-McLeod solution to the Painleve II equation 

q" = 2q 3 + sq, (A.5) 

which is singled out by the condition q(s) < for all seR. The Hastings- 
McLeod solution has the asymptotics q(s) = — Ai(s) for s —>■ oo and u(s) = 
— (— s/2) 1 / 2 for s — > — oo. (|A.3J) and f)A.4|) have to be solved with the initial 
condition 

a(s, 0) = -b{s, 0) = exp ( J ds'g(s'))- (A.6) 
The Baik-Rains scaling function is defined through 

5 , BR,(s,tf)= / ds'a(s' ,w)a(s' , —w). (A. 7) 

Proposition A.l. With the above definitions 

g(s,w) = g BR (s,w). (A.8) 

Proof. We fix w > 0. We will establish that 

d 

—g(s,w) = a(s,w)a(s,-w). (A.9) 

Then 

g(s,w) = g BR (s,w) + c. (A.10) 

Now, by construction, 

F w (s) and F* R (s) = ^-(F GlJE (s + w 2 )g BR (s + w\w)) (A.ll) 
are distribution functions with mean zero. From (jA.lOJ) we infer 

/ sdF w (s) = [ sdF* R (s)+c[ dss^F GVE (w 2 + s). (A.12) 

Since Fgue(s) is also a distribution function and ^Fgue(s) has a fast decay 
at infinity, (jA.12|) amounts to = — c and thus c = 0. 

To check (jA.9|) we will differentiate g = e w ^ 3 g. It is convenient to follow 
the scheme in j^J, Section 4.2, according to which 

/oo /*oo 
dx J dy Ai(x)p s (x,y)ty± w (y). (A.13) 
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Remark that for w > one can rewrite *ff w (x) = J° dye wy Ai(x + y) using 
We differentiate as 

— g(s,w) = e ws ^ w (s) - $ s (s) J dy Ps (s,y)^ w (y) 

POO 

-V w (s) / dx$ s (a;)p s (x, s) 



,00 ,00 9$ / V 

+ y da; y dy ^ J Ps{x,y)^M 

+ J°° dx J™ dy$ s (x) dpsi *J y) y w (y)- (A.14) 

The central identity for the proof is, see |H| I24j . 

— p s (x, y) = -— y) - y) - Q(x)Q(y) (A.15) 

as ax ay 

with 

/■oo 

Q(z) = / dyp s (x,y)Ai(y). (A.16) 



We insert flA.15|) in (jA.14j) and integrate by parts to eliminate the terms 
T^p s (x,y) and -^p s (x,y), which requires the derivatives 

d ~ 

■k-^ w (x) = M{x)-w* w (x), (A.17) 
ox 

— $ 8 (x) = w$ s (x) - Ai(x) J dykx{y)e wy - — § s {x). 
In the end only four terms remain, which can be assembled as 



■j^g(s,w)=(l-J dx^ w (x)Q(x)^J dxAi(x)e wx +£ dx$ s (x)Q(x)) . 



(A.18) 

The first factor in (jA.18|) equals a(s,w). To prove Proposition IA.1I it 

1 3 

thus remains to establish that the second factor equals es w a(s,—w). Let 
ip- w (x) = e wx . Then 

$ s = K M i/>-v, + (1 - K Ai )P s ^ w - P s ^_ w (A. 19) 

and 

Q(x) = (1 - PsKmPs^Ps Ai(x). (A.20) 
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Thus the second factor in (jA.18|) writes 

Qs Ai) + ($ a , P S Q) = faK,, Q s Ai) + (K A ^ W , P S Q) - (0_ w , p s q) 

+(V>-™,p s (i - x Ai )p s (i - p s k m p s )- 1 p s Ai) 

= (V>_ tt ,Ai) + (tf Ai ^_ w ,P,Q) 

-(if>- v ,,P.Q), (A.21) 

since the last term in the middle part equals ('0— u>, P s Ai). According to 

(jEH 

e^ 3 a(s, -w) = e^ 3 - (e^ 3 ^_„,, P,Q) (A.22) 
and by ()D.3j) it follows that e^™ 3 = Ai). Moreover, 

POO 

e^^^x) = ~ ^ / dye- wy Ai(x + y). (A.23) 



An explicit computation shows that e^^ J °° dye wy Ai(x + y) 
(KaiV'-uJOe)- Thus 



e^ 3 ^_ w = ^-^ - K M ij. w . (A.24) 
Inserting (IA.24I1 in (IA.22D one establishes that (IA.21jl equals e^ 3 a(s, 

□ 



B Determinantal fields: proof of Proposi- 
tions EE21 and HO 

B.l No boundary sources 

To prove Proposition 13.31 we find it computationally convenient to approx- 
imate the exponential distribution through a geometric one. Then w(i,j), 
i,j > 1, are independent random variables with P({u> (i, j) = n}) = (l — q)q n , 
< q < 1. The RSK construction, as explained in the main text, can 
be carried through with minor modifications, compare with [9 . The lines 
j i— > hi(j, t) take values in Z and are pinned as hi(±T, r) = £, £ = 0, —1, . . .. 
The weight of a jump of size 5 is (^/q)^- Let us denote the corresponding 
point random field by n), \j\ < t, n 6 Z, i.e., 

I 1 if there is an £ such that ht(j, r) = n, * . 

<j>%(j,n) = < (B.l) 
otherwise. 



42 



It is determinantal and, at equal times j = 2d, r = 2m, 

M 



E(n#(2d,n fc )) = det {KA n ^ n ^) i<k,k'< M - ( R2 ) 



fc=i 



We will first compute the kernel d and then show that in the exponential 
limit, lattice spacing e, q = 1 — e, 

lime-'K^de-'x], L^Vj) = K m>d (x,x') (B.3) 

for x, x' > 0. 

We have to compute d , for which we use the Fermion formalism of [2*U] . 
Let [— N, —N + 1, . . . , N] = A^v and T be the Fermionic Fock space over 
^(Aat). If A is a linear operator on £ 2 (A N ), then T(A) denotes its second 
quantization as an operator on T . Let PH be the projection onto [-N, . . . , 0] 
and let Qn be the corresponding Slater determinant. We set, as operators in 

I U tor % < j, 

and Tl = (Tl)*. T% N , Tl' N is the restriction of T|, Tl to £ 2 (Atv). Finally 
let a(j) be the Fermion field with index j, \j\ < N. Then 

K« md (i,j) = lim ^(^,r(T^T^) m+ V( J >wr(T^ iV ) m -^)^ 

(B.5) 

with the normalization 

Z N = (Q N , T(Tl' N Tl> N ) 2m n N )r. (B.6) 

Here (•, denotes the inner product in Fock space. Working out the limit, 
see EDI, results in 



K q m4 (i,j)= ^l) m+ \T q )- im - d) P^(Tl)- {m+d \Tl) m - d )(i,j), (B.7) 

where P_ = limjv^oo P^ projects onto Z_ = (...,—1,0]. 

In Fourier space T+ is multiplication by (1 — y /qe~ lh )~ 1 = Tl(k) and Tl 
multiplication by (1 — ^fqe lk )^ 1 = Tl(k). The rescaling in (|B.3J) amounts to 
replacing q by 1 — e and k by ek. Then 

limefl- £ {ek) = (1 + \ik)-\ \\mef^T e (ek) = (1 - \ik)~ x . (B.8) 



By inserting the limit (|B.8|) in (|B.7)1 . the claim ()B.3j) follows 
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B.2 Boundary sources 

We add boundary sources through the random variables w(0, 0), w(J, 0), 
w {Q,j), 3 > 1- They are independent and geometrically distributed according 
to 

P({ w (0,0) = n}) = (1 -a(3)(a(3) n } 
P({w(j,0) = n}) = (l-a^)(a^) n , 

P(Wo,j) = n}) = (l-zV^X/WN (b.9) 

with < a/5, a^g, < 1. 

The corresponding random field 0^ Q( g is again determinantal with defining 

kernel i^ 2 ^ , where we set r = 2m + 1 in accordance with the convention 

of Sectional K^™p at equal Fermionic time 2d is computed by the same 
method as in Section IB.ll In particular we first restrict the height lines to 
the interval [— N, . . . , N] and then take the limit N — > oo. Let 

f a (j) = a j , a(f a ) = J2fa(j)a(j) (B.10) 

jew 

and let be the ground state vector with sites — N, . . . , — 1 filled and sites 
0, . . . , N empty. Then a*(f a )Q]^ gives the correct weight to the jump at the 
right boundary, correspondingly for the left boundary. Therefore the defining 
kernel is given through 

K%™ +1 \2d,i;2d,j) = lim J-<£^, a(f p )T(Tl> N T^ N ) m+d (B.ll) 

a*(z)a(j)T(T^ N T^ N r- d a*(f a )Q~ N ) T 

with the normalization 

Z N = (n~, a{f p )T{T^ N T^ N Y m a*{f a )n- N )?. (B.12) 



We note that 
and 



r{A)a*(f)=a*(Af)r{A) (B.13) 

Tlf a = (1 - a^yfqY 1 }'^ y/q < a, 

Tifp = (l-/V5)~ l //», /?< * (B.14) 

v 
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Therefore 

K^ +1 \2d,i;2d,j) = Km J—(Q-,r(Tl' N T«: N r+ d (B.15) 

^ N,ap 

a{f p )a*{i)a{j)a\f a )T{Tl> N TfT- d n-) r 

with 

Zy >0/J = (f^, r(r' JV T^ iV ) m+d a(^)a*(/ Q )r(r' Ar T^ iV ) m -^)^. (B.16) 
Following the steps of Section IB.1| one concludes that 

lim Z N , a/3 = Z aP = (a(3)- l (f a , (1 - K^Jfp)^ (B.17) 

At this stage it is of use to recall a general property of quasifree states 
on CAR-algebras. Let A be a CAR algebra indexed by Z and let u be a 
quasifree linear functional (a state) on A, uniquely defined through 

u(a(j)) = co(a*(j)) = 0, ou(a*(j)a(j)) = K(i,j) (B.18) 

with K a positive linear operator on £ 2 , \\K\\ < 1. Let f E £2 an d a(/) = 
Ejez; fU) a U)- For /<> /r e ^2 we define 

Z = oo(a(f e )a*(f r )) = (f r , (1 - K)f t ) ta < 00 (B.19) 
and a linear functional u through 

Co(A) = ^u(a(f e )Aa*(f r )), A e A. (B.20) 

Then 0(1) = 1 and is again quasifree with 0(a(j)) = 0(a*(j)) = and 
covariance 

0(a*(i)a(j)) = K(i,j) + 1(1 - *T)/,(z)(l - tf)7r0')- (B.21) 

With the results from Section IB.1| the limit state (|B.15jl is a quasifree 
linear functional precisely of the form ()B.20j) with 

see jHH), and f e (j) = ft, f r (j) = a\ ^/q<a,f3< ^q~\ 
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We take the exponential limit by choosing 

1 1 

q—l — e, a = l — ea, /3 = 1 — eb, —-<a,b<-, a + b>0. (B.23) 
Then 

lim h-eadx/ej ) = Va(z) = e~ ax (B.24) 

and, for x, y > 0, 

^K^me-^YMAe-^) (B.25) 
= K m>d (x, y) + — (1 - ^m,d)^6(»)(l - K m4 )*i) a (y) 

with Z a; 6 given in (|H.24j) . This completes the proof of Proposition IH.HI 
B.3 Proof of Proposition 15^1 

Proof of Proposition VJ. 6 A By Proposition 10. II K myd is a similarity transform 
of the Laguerre kernel. Therefore K m4 = (K m4 ) 2 , ||iT m ,<z|| = 1, P u K m)d P u is 
trace class, and all eigenvalues of P u K myd P u are in the interval [0, 1]. Thus 
we only prove that, if it > 0, 1 is not in the spectrum of P u K mtd P u , which 
is accomplished by reductio ad absurdum. Assume that ip G L 2 (R + ) is an 
eigenfunction for eigenvalue 1, 

P u K m>d P u iP = (B.26) 

Then ip(x) = for x G [0, it). On the other hand, 

IMI < \\K m , d P u ^\\ < (B.27) 

Hence H-K'm^-Pu'*/'!! = IIV'II an d, since H-Km^-full = 1, one concludes that 

K mid P u i> = tfj. (B.28) 

Therefore ip is of the form (finite polynomial) x e _X//2 , which cannot vanish 
identically on [0,u). Thus the contradiction. 

To establish the second claim we use that ip a is eigenfunction of R to 
obtain 

{P u {l-K m>d )^ a ){x) = Q{x-u)(e~ ax -Z{-a) [ dyL(x, y)^ a {y)) . (B.29) 

v Jtr- j 
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Moreover, it is easy to see that \I m)C i(z)\ < 2nC mt de~ l3z for any < (3 < 1 — p, 
C mt d being a constant (take as path Ti_ p the circle centered in 1 — p of radius 
1 — p — (3). Let us choose any (3 G (1/2 — p + a, 1 — p). Then, for u > 0, 

\{P U {1-K m4 )ij a ){x)\ < Q(x)e~ ax + (B.30) 

+&(x)Z(-a)C m4 e-^-^ +p)x [ dye^ +p -^- a) . 

(3 — ~ + p > 0, because a > 0, and /3 + p — ^ — a>0, because (3 > ~ — p + a. 
Thus P u (l — Km^ipa £ £ 2 (1R) with norm uniformly bounded in w. The 
second part of ()3.22j) is treated similarly. ip a is eigenfunction of L*, thus 

(P u (l-K m , d )*i> a )(x) = 6(x- M )(e- aa; -Z(a)- 1 ! dyR(y, x)e~ a " > j . (B.31) 

J IR_ 

l-^m^C^)! < ^C m ^e~^ z for any < (3 < p. Thus, by choosing (3 G (p + a — 
1/2, p), for it > 0^ 

|(P U (1 - K m , d )Va)(^)| < &(x)e~ ax + (B.32) 

+0(x)Z(-a)(7 m , d e-( /3+ 5-^ / dye^*^, 

JlR- 

which, by the choice of f3 and since a > 0, implies P u (t — K m ^)*ip a G L 2 (R) 
with norm uniformly bounded in u. 

□ 

C The Laguerre kernel 

Let L„ be the standard n-th Laguerre polynomial of integer order a, a > 
|16j . Chapter 5.3. The Laguerre polynomials are orthogonal on R + relative 
to the weight x a e~ x as 

L dxx " e " ((i^of) Li ° >(:c) ((^W) ^ w = (ai) 

The Laguerre kernel is the orthogonal projection onto the first n Laguerre 
polynomials and is given by 

n— 1 .| 

Ki a \x,y) = T-j^Lf^Lf^x^y^e-^e-y^ (C.2) 
(j +a)\ 

for x, y > 0. 
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Proposition C.l. With definitions and \C. 2i) one has 

K m _ d (x, y) = K mid (x, y) (C.3) 

and 

K m4 (x,y) = K^(x,y)Q d (C.4) 

for < d < m and x, y > 0. 
Proof. Let us define 

m = {\-^k) m - d {\+lk)- {m+d \ 

g r (k) = (\-ik) m+d (\ + ik)- {m - d) . (C.5) 
We set a = 2d, n = m - d - 1. By [IE], p. 244, it follows 



dke**g t (k) = gt(x) = (-1)" ^ .. e^l - d x )x a L^\x)Q{x) (C.6) 
r (n + a)! 



and 



/ dke- ikx g r (k)=g r (-x), (C.7) 



with 

Orix) = (| + - ^) d+1 (-l)" l+1 e-^ 2 L^ 1 (x)e(x). (C.8) 

^ is supported on [0, oo) with a discontinuity at x = 0. As a distribution g r 
is supported in [0, oo) with the singular part concentrated at {x = 0}. 
By [IE], Section 5.5.2, one has the identities 

(l-dM a) = ~dA% &Lto = -L%?>. (C.9) 
Using them repeatedly in (|C.6|) and ()C.8|) yields, for x > 0, 

= (-lT +1 T J ^e-*/ 2 (x a d x LW 1 (x) + ax^L^), (CIO) 

and 

O-On) = (-l^e-^^Ca:). (C.ll) 
With these notations the integral kernel from ()3.13|) is expressed as 

^ T m+d T -(m-d) p _ T -(m+d) T rn-dj ^ ^ ^ ^ 

dwg e (x - w)g r (-(w -y)) = / dwg e (x + w)g r (y + w). 

LR Jr + 
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Note that for x, y > only the regular part of g r is used. We insert (IC.lOjl 
and (IUXT| in (RTHJ) . Then (ITIIJ) of Proposition RTT1 amounts to 



x ((x + w) Q 9„,4 Q i(x + w) + a(x + w) a - l L^\y + w)). (C.13) 

We check recursively by setting the left hand side as 52" =0 iVO + a V-Bj 
and the right side as n\/{n + a)\A n . Then (IC.13J) is equivalent to 

n\ (n — 1)! n! 

Ao = 5 °' (^)! n " (n - 1 + a)! = (^)! n = 1 ' 2 ' — 

(C.14) 

A = -Bq amounts to a partial integration. For the second equality we write 



A n - ^—^-A n _i = / dwe- w (x + w) a (L' n+1 (x + w)L' n+1 {y + w) 

n J R+ V 

_!l±^ L ' n ( x + + w ))+ [ dwe- w a(x + w)*- 1 

x (^L n (x + tu)Z/ n+1 (y + u>) I, n _i(ar + iw)IrJ l (3/ + iw)J, (C.15) 

omitting the superscript a. In the second integral we use the identities 
L' n+1 = L' n — L n and xL' n = nL n — {n + a)L n _i. Then the terms combine as 

A n ~^-^-A n _i = I dwe~ w (x + w) a [L n {x + w)L n (y + w) 

n JTR+ 

-L' n (x + w)L n (y + w)- L n (x + w)L' n (y + w)) 

— I dwe~ w a(x + w) a ~ 1 L n (x + w)L n (y + w) 
Jr+ 

= - / dw— (e~ w {x + w) a L n {x + w)L n {y + w)) 
Jn+ dw y 

= B n , (C.16) 
which is the recursion relation (jC14|) . □ 

D Some useful relations 

D.l Two representations of the Airy functions 

For any cr < 0, define the path 7^ = 17 + zR. Then 

Ai(z) = -L / d£e-^ 3+ *« (D.l) 
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and for any /x > 0, with 7^ = /i + iH, 



M(z) = — I d^ 3 -^. (D.2) 



2iri 



7m 



One can deform 7^ so that it goes from ooe~ 2m ^ to ooe 2m ^ and crosses the 
real axis at a. In this case the paths will be denoted by 7*. Similarly the 
deformation of 7^ goes from ooe _Mr / 3 to ooe t7T ^ 3 and is denoted by r y r . 
A formula which will be employed later is 

M(w) = [ dye wy Pd((3 + y)=e^- pw (D.3) 



K 



valid for all u> > (for w = as improper Riemann integral). To prove it 
one derives the differential equation 

^± = M(w)(w 2 - f3) (D.4) 
aw 

by integrating twice by parts. ()D.3j) follows from (|D.4)1 and the initial value 
M(w = 0) = 1. 

D.2 Two integrals around a pole 

Let be the path enclosing z = w and anti-clockwise oriented. Then 

/ d £ ^"7 = Res — ;£ = w ]=-8 

(D.5) 

and 

s. 



f e -^+^w i +(w 2 +s)((~w) / e -^ 3 + ^w s + (w 2 +s)(Z-w) \ 

X. df m - ^ = R H — <r*? — ;f = w ) 



D.3 Two equivalent expressions for S v 



(D.6) 



Let us consider w > 0. The representation of the Airy function (jD.2|) allows 
one to write 

POO POO 

S w>s = s + e ~\ w3 dx dy Ai(w 2 + s + x + y)e w{w2+s+x+y) 
Jo Jo 

poo POO -I p 

= s + e^ w3 / dx / dy— / d £e^V w2+s+aH *-« (D.7) 
Jo Jo 27 ™ J% 
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By choosing /i = 2w, one obtains an integrand which is absolutely inte- 
grable. Thus by Fubini's Theorem we are allowed to exchange the order of 
the integrals and compute first the one in (x, y) with the result 



Let r„, be as in (ID. 51). then 



Sw ' S = 2m(w-0> + J r ^ 2m(w-0 2 

e §? 3 -f™ 3 +(™ 2 +s)(>-£) 

d £ V~( 7V2 • ( D - 9 ) 

2m (w — f Y 

On the other hand, 
/ dxdy$ W!S (x + y) = [ dxdy— [ d^ 3 -^ 3 e ^-0(^+s+x+y) 

d £ o-„V„ . ^2 = ( D ' 10 ) 

w/2+m. 



2m{w - 2 



by deforming the integral on w/2 + iH to 7^ 



w/2- 



References 

[1] J. Baik, G. Ben Arous, and S. Peche, Phase transition of the largest 
eigenvalue for non-null complex sample covariance matrices, Ann. 
Probab. 33 (2005), 1643-1697. 

[2] J. Baik and E.M. Rains, Limiting distributions for a polynuclear growth 
model with external sources, J. Stat. Phys. 100 (2000), 523-542. 

[3] J. Baik and E.M. Rains, Symmetrized random permutations, Random 
Matrix Models and Their Applications, vol. 40, Cambridge University 
Press, 2001, pp. 1-19. 

[4] H. van Beijeren, R. Kutner, and H. Spohn, Excess noise for driven 
diffusive systems, Phys. Rev. Lett. 54 (1985), 2026-2029. 

[5] F. Colaiori and M.A. Moore, Numerical solution of the mode-coupling 
equations for the Kardar-Parisi-Zhang equation in one dimension, Phys. 
Rev. E 65 (2002), 017105. 



51 



[6] P.L. Ferrari, Polynuclear growth on a flat substrate and edge scaling of 
GOE eigenvalues, Comm. Math. Phys. 252 (2004), 77-109. 

[7] P.L. Ferrari and H. Spohn, A determinantal formula for the GOE Tracy- 
Widom distribution, J. Phys. A 38 (2005), L557-L561. 

[8] D. Forster, D.R. Nelson, and M.J. Stephen, Large- distance and long-time 
properties of a randomly stirred fluid, Phys. Rev. A 16 (1977), 732-749. 

[9] T. Imamura and T. Sasamoto, Fluctuations of the one- dimensional 
polynuclear growth model with external sources, Nucl. Phys. B 699 
(2004), 503-544. 

[10] K. Johansson, Shape fluctuations and random matrices, Comm. Math. 
Phys. 209 (2000), 437-476. 

[11] K. Johansson, Discrete polynuclear growth and determinantal processes, 
Comm. Math. Phys. 242 (2003), 277-329. 

[12] K. Kardar, G. Parisi, and Y.Z. Zhang, Dynamic scaling of growing in- 
terfaces, Phys. Rev. Lett. 56 (1986), 889-892. 

[13] J. Krug, P. Meakin, and T. Halpin-Healy, Amplitude universality for 
driven interfaces and directed polymers in random media, Phys. Rev. A 
45 (1992), 638-653. 

[14] T.M. Liggett, Coupling the simple exclusion process, Ann. Probab. 4 
(1976), 339-356. 

[15] T.M. Liggett, Stochastic interacting systems: contact, voter and exclu- 
sion processes, Springer Verlag, Berlin, 1985. 

[16] W. Magnus, F. Oberhettinger, and R.P. Soni, Formulas and theorems 
for the special functions of mathematical physics, Grundlehren Band 52, 
Springer Verlag, Berlin, 1966. 

[17] A. Okounkov, Infinite wedge and random partitions, Selecta Math. 7 
(2001), 57-81. 

[18] M. Prahofer, Stochastic surface growth, Ph.D. 

thesis, Ludwig-Maximilians-Universitat, Miinchen, 

http://edoc.ub.uni-muenchen.de/archive/00001381 , 2003. 

[19] M. Prahofer and H. Spohn, Current fluctuations for the totally asymmet- 
ric simple exclusion process, In and out of equilibrium (V. Sidoravicius, 
ed.), Progress in Probability, Birkhauser, 2002. 



52 



[20] M. Prahofer and H. Spohn, Scale invariance of the PNG droplet and the 
Airy process, J. Stat. Phys. 108 (2002), 1071-1106. 

[21] M. Prahofer and H. Spohn, Exact scaling function for one- dimensional 
stationary KPZ growth, J. Stat. Phys. 115 (2004), 255-279. 

[22] T. Sasamoto, Spatial correlations of the ID KPZ surface on a flat sub- 
strate, J. Phys. A 38 (2005), L549-L556. 

[23] H. Spohn, Excess noise for a lattice gas model of a resistor, Z. Phys. B 
57 (1984), 255-261. 

[24] C.A. Tracy and H. Widom, Level-spacing distributions and the Airy 
kernel, Comm. Math. Phys. 159 (1994), 151-174. 



53 



